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It has recently been proposed, with the help of numerical investigations, that fluorescence cor- 
CNj , relation spectroscopy at variable observation area can reveal the existence of a meshgrid of semi- 

permeable barriers hindering the two-dimensional diffusion of tagged particles, such as plasmic mem- 
^ ' brane constituents. We present a complete theory confirming and accounting for these findings. It 

, enables a reliable, quantitative exploitation of experimental data from which the sub-wavelength 

■ mesh size can be extracted. Time scales at which fluorescence correlation spectroscopy must be 

' performed experimentally are discussed in detail. 

<n : 

r ~~ 1 ' I. INTRODUCTION 

o ; 

The rapid development of experimental techniques for the exploration of the structure of cellular membranes at the 
sub-micrometric scale emphasizes the need for suitable tools to interpret the large amount of data collected. Among 
the many techniques routinely used in biology or biophysics laboratories, Fluorescence Correlation Microscopy (FCS) 
0, Q is very promising because it has excellent temporal dynamics: the auto-correlation function can be built for 
time-scales going from nanoseconds to minutes. Even though its spatial resolution is a priori limited by optical 
diffraction (typically 200 nm), recent physical improvements have broken this Rayleigh limit [!, Q. Their spatial 
resolution can go down to a few tens of nanometers, thus making accessible the characterization of nano-domains in 
£h cell plasma membranes Q . 

With this powerful tool, it is now crucial to be able to identify, in the experimental data dealing with real biological 
systems, the signatures of the different diffusion modes @, @, 0|- Indeed, molecules are in constant movement in live 
systems because the cohesion energies between biological molecules are generally of the same order of magnitude as 
the thermal energy fc^T. Thermal agitation is omnipresent at this scale and must be taken into account in the analysis 
of biological processes. In particular, thermally induced diffusion certainly plays an important role in many biological 
functions at the molecular scale, and there is everyday new evidence that characterizing these functions cannot go 
without the understanding of the subtle dynamical organization of cell constituents. As far as proteins or lipids 
embedded in the cellular plasma membrane are concerned, their diffusion is furthermore affected by heterogeneities of 
| the membrane or of its close environment. For instance, they can be confined or partially confined in micro-domains. 
\Q . This confinement is supposed to enhance the interactions with co-confined constituents and plays a role in the biological 
functions of membrane molecules. The processes responsible for confinement have been widely investigated in the last 
. decade and several mechanisms have been proposed to account for it, depending on the biological systems (§. Itj ITol. [Tl| : 
y—i ' see also the review Without giving greater importance to any one of these mechanisms, we focus here on the 
f — case where diffusion is hindered by linear semi-permeable obstacles ("fences"), whatever their nature, constituting a 
meshgrid with a typical mesh size (or mesh parameter). Cortical cytoskeleton filaments have been demonstrated to 
play such a role by slowing down the long-term diffusion of membrane constituents [H, [l3j , with a mesh size of several 
tens or few hundreds of nanometers. The diffusion appears to be confined at short times and the ensuing long-term, 
macroscopic diffusion coefficient, Dm, is much smaller than the short-term, microscopic one, D^. 

In this context and following previous works on Fluorescence Recovery After Photo-bleaching (FRAP) at variable 
radius [3, EH Ho ]. P.-F. Lenne and his collaborators have recently demonstrated that the use of FCS at variable 
observation area enables one to distinguish between several diffusion modes d H, 0). They proved, with the help 
of numerical simulations, that the three following diffusional behaviors have a different signature: free diffusion; 
diffusion on a membrane containing disjoint sub-micrometric domains where membrane molecules are temporarily 
trapped and diffuse more slowly, with the idea of identifying raft-like structures; diffusion on a sub-micrometric 
meshgrid constituted of joint domains, inside which the molecules exhibit confined diffusion, but that are separated 
by semi-permeable barriers over which the molecule can jump with a certain probability only, as discussed above. 
Our goal in the present paper is precisely to propose a rigorous analytical calculation of the two-time fluorescence 
autocorrelation function in the case of molecules (proteins or lipids) diffusing on such a meshgrid, thus granting the 
request by M.J. Saxton: "The theory remains to be done" [7|. We shall prove that if ii/2 is the time at which the 
FCS correlation function gets reduced by a factor 2, if w is the waist (defined below) of the Gaussian observation 
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profile, then 



t 1/2 ^——(w 2 -a 2 /12), (I) 



where a is the mesh size. Experimentally, the measure of ti/2 at several waists enables one to infer the value of a. 



The two-dimensional observation profile used in FCS experiments is usually close to Gaussian [17j, whether it is a 
focused laser beam or the two-dimensional section of a confocal microscope observation volume. For these reasons, we 
shall compute the theoretical correlation function in this Gaussian context. However, because we think that a square, 
uniform profile is better adapted to introduce the basic concepts useful in the calculation, we will first focus on this 
academic uniform case before tackling the Gaussian one. In recent works, nanometric apertures enabling the use of 
observation profiles much smaller than the Rayleigh limit have been developed 0, Q ■ A priori, these observation 
profiles are closer to a uniform disc than to a Gaussian, even though the subtle optics through a nano-aperture smaller 
than the light wavelength deserves a more thorough analysis (see the discussion in Ref. Q). However, as it will be 
discussed below, the technique developed in the present paper does not allow the analysis of a uniform, circular, 
observation profile. We shall restrict it to square and Gaussian ones, and we shall see that our conclusions will be 
qualitatively equivalent in both cases. 

After revisiting in the following Section [IT] the standard background, among which free diffusion, the central Sec- 
tion IIIII will be devoted to the main results of this paper concerning the diffusion on a meshgrid. Most numerical 
results obtained in previous references will be explained, including the short-time shape of the correlation functions. 
Finally, the last section ITVl contains concluding remarks as well as a discussion with an experimental perspective. The 
two appendices contain technical material useful in Section IIIII 



II. BACKGROUND, DEFINITIONS AND FREE DIFFUSION 



Here we briefly recall some generalities about FCS. As justified in the Introduction, the context is the diffusion 
of independent particles in a quasi- two-dimensional medium (e.g. proteins or lipids in a lipidic bilayer), tagged by 
fluorescent labels that are observed thanks to a focused light source illuminating a small region of the sample, the 
area of which varies typically between a fraction of square micrometer to several ones. 

In particular, we explain here how to compute the correlation function for freely diffusing particles. This simple 
calculation will be useful to understand more complex ones in the following sections. Furthermore, we must recover 
these free-diffusion laws for diffusion on a meshgrid at large time or space scales, which will be a test of the validity 
of our calculations. As discussed previously, we first focus on a square observation profile before considering the 
more physical Gaussian case. We also show how the two-dimensional problem can be written as a product of two 
one-dimensional ones, thus reducing our investigations to ID diffusion. 



A. Correlation function 



In FCS experiments, the two-time auto-correlation function, g^(t), indicates the average degree of correlation 
between measures of the number of photons, n(s) and n(s + t) collected by the detector, at different times s and 
s + t P, 0. The decorrelation depends on the diffusion coefficient (s) of the tagged particles. Each measure of n is a 
number of photons detected per sampling time At. More precisely, the correlation function is defined as follows: 

(n(s)) 2 

The brackets mean an average over a large amount of realizations of the system (ensemble average). Because of 
translational symmetry in s at equilibrium, the previous quantity is independent of s. However, owing to the ergodic 
theorem, it can also be seen as an average over successive times s in a same, very long trajectory (temporal average). 
When needed, we shall adopt the latter point of view in the following. 

To detect fluorophores in the observed area, an excitation light is used, with a profile I(r) = I(x, y) depending on 
the two-dimensional position r = (x, y) on the sample. Then the number of collected photons at time s is 

n(s) = At J dxdy I(x,y)Qp(x,y, s), (3) 

where the numerical factor Q is the product of the absorption cross section of one fluorophore by the fluorescence 
quantum yield and the fluorescence efficiency 0, H3] and p(x, y, s) is the number of fluorophores per unit area at time 



3 



s, whether they are observed or not during At: p depends on sample and realization only, and not on the illumination 
profile. A straightforward calculation [2| shows that at the thermodynamic limit, g^ 2 \t) has the typical form: 

9 {2] {t) = l + -± T -C 1(r) {t). (4) 

In this expression, (p) is the ensemble average of the density of fluorophores in the system. From now, we assume 
without loss of generality that (p) = 1. Furthermore, S e s is an effective observation area: 

s ^ _ [J dxdy I(x, y)} 2 
j dxdy [I(x,y)} 2 

As for Cj( r )(£), it is the core of the auto-correlation function. With our choice of normalization by S e s(p), it decreases 
from 1 at t — to at large times. Its precise form depends on the observation profile I(r). The prefactor l/(S c g(p)) 
is equal to the usual one, l/N, found in the literature @, H, H, E3- 

At the experimental level, the measure of g^ (t) on time-scales ranging from nanoseconds to minutes provides a 
very good estimate of C7( r ) (t). An important quantity that is extracted from FCS data is the diffusion or decorrelation 
time, ti/2, defined by 

C I{r) {t 1/2 ) = \. (6) 

By a scaling argument, one can already foresee that t\i 2 will be proportional to the ratio of the observation area and 
a diffusion constant. The aim of the present work is to precise this rough argument and to explicit what additional 
information can be extracted from FCS experiments. 



B. Normalized observation profiles 



As justified in the introduction, we study in this paper two different two-dimensional illumination intensity profiles 
(observation profiles) I(x,y): a square, uniform one and a Gaussian one. We normalize them by J I(x,y) dxdy = 1, 
without loss of generality since the normalization disappears in Eq. The square profile is 

I( X >y) = j2M~L/2,L/2]{x)l[-L/2,L/2]{y), (?) 

where 1[_l/2,l/2] is the characteristic function of the segment [— L/2,L/2], equal to 1 on the segment and to 
anywhere else. One calculates that S e s = L 2 . As for the normalized Gaussian profile, 

. , 2 / 2(x 2 + y 2 )\ 

I(x,y ) = — exp -- i ^ LL ), (8) 

irw \ w J 

where w is the profile waist (twice its mean standard deviation (a; 2 ) 1 / 2 ). Furthermore S e f[ = nw 2 . 



C. From two to one dimension 



In this paper, we focus on diffusion in two dimensions. But it appears that the problem can be reduced to a 
one-dimensional one because: (i) the square meshgrid is a product of two one-dimensional grids of mesh parameter 
a; (ii) the diffusion operator - the Laplacian in the free-diffusion case, or a more complex Fokker-Planck one fljl ] in a 
meshgrid, see below - is a tensor product of two one-dimensional diffusion operators; and (iii) the observation profile 
can be written as a product of two one-dimensional ones, both in the square and in the Gaussian cases: 

I{x,y) = h(x)I 2 {y), (9) 

where I\{x) and ^(y) ar e normalized. If the observation profile is uniform on a square, then I\{x) and I 2 {y) are 
proportional to the characteristic function of a segment of length L. Such I\(x) or I 2 (y) will be referred to as "uniform 
profiles" in the following. If the two-dimensional profile is Gaussian, then Ii(x) and I 2 (y) are also Gaussian, with the 
same waist as I(x,y). To simplify, I\{x) will also be written I(x) when this notation is not ambiguous. 
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In these conditions, it is proven |17( that the correlation function Cjm above is nothing but the product of two 
one-dimensional correlation functions C(i), also depending on L or w. The initial problem reduces to the study of 
one-dimensional diffusion on a lattice of equally spaced punctual barriers that can be passed with a certain probability. 
In the cases under interest here, the observation profile is either uniform of length L or Gaussian of waist w. The 
only difference is that the decorrelation time t 1 / 2 , a function of L or w, is now the solution of 

C(h /2 ) = -J=. (10) 

We must also take care of tracking correctly the prefactor 1/S e s coming from our choice of normalization of C/( r ) (the 
case of the factor l/(p) has been fixed by setting it to 1): a prefactor l/^/S c s is distributed to each one-dimensional 
C (see below). 

Note that we could factorize in the same way a two-dimensional problem where the grid has a two different mesh 
parameters a x and a y ; or where the observation profile has two different characteristic lengths in the x and y directions 
(L x and L y or w x and w y ). For simplicity's sake, we shall focus on the "isotropic" case. 

By contrast, a circular observation profile cannot be written as a product of one-dimensional profiles. It is for this 
reason that this case cannot be addressed by using the tools proposed in the present paper. Note also that the previous 
formal decomposition in a product of one-dimensional problems can be interesting from a numerical perspective @ 
since it reduces the computational complexity. 



D. Calculation of the free-diffusion correlation function 



Here we give a rapid summary of calculations that can be found in former references dealing with FCS, e.g. 
Refs. [HQ- In this special case, the function C, the diffusion coefficient and the decorrelation time t 1 / 2 are denoted 
respectively by T, Q and to (thus T(<o) = l/\/2)- We also introduce the adimensional variable x: 

2v^t , 2Vtti. 

x = — - — (resp. x — ) (11) 

L w 

for a uniform (resp. Gaussian) observation profile. Since, by scaling arguments, to will be of order L 2 /fi or w 2 /fi, 
the value of x corresponding to to (denoted by xq) will be of order 1 in time intervals of interest. 

If I(q) = 1/V27T J cxp(—iqx)I(x)dq is the Fourier transform of I(x), then the correlation function T(t) depends on 
the observation profile as follows [17] : 

roc 

-Q.q 2 t 



r(t) = V^eff/ dq\I{q)\'e- Uqt . (12) 

J — oc 

In this equation, the only contribution of the underlying free diffusion is contained in f2. Note that T(0) = 1, 
as required, by definition of S e g and because I is normalized. For a uniform observation profile, one gets after 
integration 

Tlx) = erffl/x) - 4=(1 " exp(-l/a; 2 )), (13) 

where erf is the usual error function [l9j . The graph of T is displayed in Figure [T] By definition, xq is defined by 

T(x ) = 1/V2. The unique solution is xq ~ 0.520, a transcendental real. At the leading order in L, t 1 / 2 = . 
This means that the ti/ 2 vs L 2 profiles are strictly linear with slope :Eq/(4£1). 

For a Gaussian observation profile, one gets the following function, the graph of which is also given in Figure [TJ 
r(.x) = l/Vl + x 2 . In this case Xq = 1, a much simpler solution, and the decorrelation time obeys ti/ 2 = L 2 /(4£l), 
and r(t) can be written 

r(t) = 1 (14) 
^ + t2 /t\ /2 

Since xq ~ 0.520 ~ 1/2 in the uniform case, L/2 and w play essentially the same role above, which means that the 
waist w can be compared to the half-width of a uniform observation profile for FCS purposes. 
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FIG. 1: Fluorescence correlations in dimension 1: Graphs of the functions T, as defined in Eqs. (|13[1 and (|14[) . for the uniform 
(upper curve) and Gaussian (lower curve) observation profiles, in log-linear coordinates. The reduced variable x is defined as 
x = 2V Dutj L in the uniform case ori = 2\fDMt/w in the Gaussian one. Both functions essentially differ by an abscissa- shift. 

III. DIFFUSION ON A MESHGRID 

Here we consider a two-dimensional system still constituted of independent particles, but the diffusion of which is 
now hindered by a meshgrid of linear obstacles. For simplicity [||, we suppose that the meshgrid is square, with mesh 
parameter a, as illustrated in Fig. [2l The edges of the meshgrid are semi-permeable walls: they generally prevent 
the motion of diffusers, resulting in free diffusion with reflecting boundary conditions in square regions. We denote 
by the associated (microscopic) diffusion coefficient. But the particles can also jump over these obstacles with a 
low probability (e.g. by thermal activation). In the following, the squares constituting the meshgrid will be referred 
to as "domains" or "boxes" , in which the diffusion is confined at intermediate times. At larger times, we expect a 
slow diffusion mode, with a (macroscopic) diffusion coefficient, denoted by Dm, and calculated below. It should be 
governed by the laws derived in the previous section, where the free diffusion coefficient fl has to be replace d by 
Dm- We suppose that time-scales are well separated, as it is usually required in such a context 0, @, [T3, El, Il6j : 
Dm <C D^. This means that the barriers are true (semi-permeable) obstacles, that have a chance to be identified as 
such experimentally. 

A question that arises at this stage is that of the position of the observation profile relatively to the meshgrid. 
For example, it has already been observed 0, @ that FCS profiles differ according to whether the center of the 
Gaussian profile falls close to a vertex or near the center of a domain. This issue will be carefully questioned in the 
following. Experimentally, it is common to perform many experiments on the same (model or live) system before 
averaging observables. This procedure is used to measure t 1 / 2 . However, the averaging process can be performed at 
two different steps, leading a priori to slightly different values of t 1 ^ 2 . Either different correlation functions C(t, uu) are 
measured before averaging them. The resulting average is denoted by C(t,w) and t 1 / 2 is defined by C(t 1 / 2 ) = 1/2. 
Or one ti/ 2 is measured for each experiment and the final value of ti/ 2 is their average over realizations. In the 
following calculations, we shall adapt the first point of view, for the sake of simplicity. However, we shall prove in 
the last section devoted to discussions that both procedure are equivalent for experimental determination of the mesh 
parameter a, the main purpose of the present work. 

A. Model, box approximation and time scales 

As it was justified above, we now focus on the one-dimensional case, that is to say the diffusion on a line with 
regularly spaced semi-permeable punctual obstacles. Anticipating on the following, we refer to Figure [3] where an 
example of correlation function C(t, L) is displayed (upper curve). It corresponds to a uniform observation profile but 
it also qualitatively matches the Gaussian case. One clearly identifies two time-scales in this graph, corresponding 
to the two successive fall-offs in the correlation function. We shall see that the first (small) fall-off, around t^, is 
associated with the "microscopic" decorrelation inside domains on the border of the observation area; The second 
(large) one, around tm, accounts for the "macroscopic" decorrelation inside the whole observation area. 
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FIG. 2: The simulated trajectory of one particle diffusing in a regular two-dimensional meshgrid of mesh parameter a — 400 nm. 
It diffuses freely in each square domain, with diffusion coefficient arbitrarily set to D M = 0.06 /im 2 s _1 , and it is generally subject 
to reflecting boundary conditions on the square edges. However, it can jump over edges with a low probability, thus resulting in 
an apparent slow diffusion at large times, with diffusion coefficient Dm = here. The circle materializes the observation 

area, in which the number of observed particles fluctuates. The correlation function gr 2 ' measures these fluctuations and their 
typical time scales. 



In order to compute the correlation function C (t) , we shall also need to separate the slow and rapid time scales in 
the averaging process. Having already identified the two typical "microscopic" and "macroscopic" time scales r M and 
tm , we introduce an additional intermediate "mesoscopic" time scale, r m , such that t m <C r m <C tm ■ When averaging 
over times t in a trajectory below, we shall distinguish between short-term averages, denoted by {-) s t and long-term 
ones, denoted by (•)/*■ This means that we first calculate averages on segments [0, r m ], [r m , 2t,„], etc. . . , then we take 
the mean value of these short-term averages. Formally, in Eq. @, (•) = ((^stjit- In other words, in the reciprocal 
space, (-} s t (resp. (-)it) means integration over frequencies larger (resp. smaller) than l/r m . 



The model used in this paper to tackle the diffusion on a meshgrid is schematized in Figure 2J Particles diffuse in 
a periodic potential U(x), of period the mesh parameter a. It is constituted of wells separated by barriers of width 
b -C a and of finite height Uq ("punctual" semi-permeable obstacles). The fact that barriers are semi-permeable 
implies that Uq ksT so that jumps over barriers, requiring an energy Uq, are rare. This argument will be refined 
below. The particle density p(x, t) is a variable indicating the probability, p(x, t) dx, of finding a particle at time t in 
the interval [x, x + dx]. By definition, at equilibrium, (p(x, £)) = (p) = 1, in conformity with our previous prescription. 

For the sake of generality, we can set two different diffusion coefficients in wells and barriers, denoted respectively 
by and The equation governing the temporal evolution of p is then the Fokker-Planck equation . 18j |. This 
equation can be solved in the present case. However, because this generalization does not bring any insight into the 
physics and because it needlessly complicates the mathematical formalism, we suppose from now that Df, = . The 
Fokker-Planck equation then amounts to the simpler Smoluchowski one With our choice of U(x), the situation 
simplifies again because a classical diffusion equation is valid inside wells or barriers, where both the potential, U, 
and the diffusion coefficient, D M , are constant. But it is not able to tackle the variations of U with x. In the special 
case under interest here, the Smoluchowski equation simplifies into the "jump condition" (l8j because there exist only 
isolated discontinuities of U: if the discontinuity occurs at x = Xq, and if AU = U (xq + s) — U(xq — s), then 

p(x +e) = exp \-^7fj P( x o - e)- (15) 

Since Uq <g; kgT, p is much smaller inside barriers than outside, as illustrated in Figure 5] 

The principle of the quasi-static box approximation scheme developed here is also exemplified in this figure. The 
equilibration (or decor relation) time-scale inside a well, denoted above by r M , is given by the resolution of the diffusion 
equation in a square well, the boundaries of which can be considered as perfectly reflecting since jumps are scarce. 
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FIG. 3: Example of fluorescence correlation function C(t, L) for a uniform observation profile of length L — 4, after integration 
on the position of the profile center relatively to the meshgrid of mesh parameter a. The time unit to is arbitrary. The diffusion 
parameters are Dm = 10 _4 a 2 rJ" 1 and = 5.10 _2 a 2 T -1 . The black curve shows the entire profile whereas the red one displays 
the long-term component only. The "microscopic" time-scale r M ~ 1 to and the "macroscopic" one tm ~ 10 5 to are clearly 
identifiable, corresponding to the two successive fall-offs in C. The "mesoscopic" time-scale r m in-between (gray vertical bar) 
separates the short-time regime on the left and the long-time one on the right. The decorrelation time tj/ a is the unique 
solution of C(ti/2, L) = \ j\f2. The knowledge of the long-term component suffices if one wants to calculate ti/%. 




FIG. 4: Quasi-static box approximation in one dimension. The boxes of width a are separated by potential barriers (in black) 
of width b <C a and height Uo ksT. The particle density p (in red) is quasi-uniform in domains because the exchange rate 
between domains is much slower than the equilibration time inside domains. The difference of filling between adjacent domains 
results in a weak density current, J = — _D M Vp, at the barrier, due to the gradient of p inside the barrier. This is the cause of 
the slow long-term diffusion. 



One gets (see Section HHD[) : 

,2 



(16) 



This means that a non-uniform distribution p inside the well will equilibrate and become uniform after a time much 
larger than r M . For the example given in Figure[3l on gets r M ~ 2.0, as anticipated. 

Now we focus on what happens inside the barriers. Since b <C a and particle densities are much smaller inside 
barriers (Eq. (|15p). the probability of finding a particle inside a barrier is very small as compared to the wells. 
Wells behave like reservoirs with regard to barriers and impose boundary conditions to p in barriers. In addition, 
the equilibration time inside barriers, r& is also given by Eq. (| 16p : Tf, = b 2 / (ir 2 D M ) . Then the diffusion equation 
becomes d 2 p/dx 2 = 1/-D„ dp/dt = (equilibrium condition), that is to say p{x) — Ax + B. Thus p is affine in 
barriers, its values on the boundaries being fixed by the filling of both adjacent boxes (see Figure [4]). The quasi- 
static approximation assumes that all time-scales in the problem are much larger than t m and r e (which is checked a 
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posteriori for self-consistency), so that the p profiles inside boxes and barriers can be considered at equilibrium. They 
can vary with time, but very slowly. 

Indeed, suppose that two adjacent boxes of indices i and i + 1 are differently filled, as in Figure |U Then a particle 
current appears through the barrier from the more filled box towards the less filled one. This current is denoted by 
Jj— n+i- It is positive if particles go from i to i + 1, and negative in the converse case. More quantitatively, the total 
particle number inside the box of index i at time t, J box /5(x,t) da;, is denoted by Wi(t). The particle density on 
the box boundary is therefore Wi(t)/a. The condition on the left-hand (resp. right-hand) boundary of the barrier 
between boxes i and i + 1 is thus p — Wi(t)/a exp(— C/o/feT) (resp. p = Wi + \(t)/a exp(— Uo/ksT)). Since p is 
affine inside the barrier, Vp = exp(— Uq/ kBT)(Wi+i(t) — Wi(t)) / (ab) and 

Ji-i+i = -D,Vp = -D,cM-U /k B T) W *+^ b W ^\ (17) 

In the following section, we shall write and solve the diffusion equation satisfied by the Wi(t). It will ensue from the 
previous expression of J. 

To conclude this section, we have demonstrated that, in the conditions under interest in this paper, the diffusion on 
a meshgrid can be treated in the quasi-static box approximation, where domain fillings, Wi(t), are uniform and vary 
slowly. With respect to the short-term and long-term averages discussed above, it means that the short-term average 
(■) s t is implicitly already realized when supposing that p(x,t) = Wi(t) in box i, since p does not display short-term 
fluctuations any-longer. Again implicitly, all remaining averages (•) are in fact long-term ones, (-)it- Finally, we 
remark that because b <C a and particle densities are small inside barriers, the probability of finding a particle inside 
a barrier is very small and will be neglected. 



B. Discrete diffusion equation and box-box two-time correlation function 

The evolution of Wi(t) ensues from the currents J^i+i and Ji_»j_i: 

>,(*) = -(J M+1 + WO - D M "■-"»- 2H '"» + W - +1 " ) - 08) 

at a z 
This is a discrete-position diffusion equation, the right-hand-side fraction being nothing but a discrete Laplacian on a 
one-dimension lattice of lattice-spacing a. Therefore, as anticipated, the long-term behavior of particles is diffusional, 
with the long-term diffusion coefficient 

a __Hs_ , , 

D M = -rD^e "b? . (19) 

The three parameters b, D^, and Uq eventually amount here to a single one, Dm- Note that the condition Dm <C D^ 
required above imposes an additional condition on b: a b 3> aexp(— U^/hsT). In practice, Uq > 5k B T suffices to 
ensure that such a value of b can exist: b ~ a/10 satisfies both conditions. 
A new box-exchange time-scale appears in this equation, 

a 2 ab _Hs_ 

re = 7r = 7T e ' BT ' ( 2 °) 
Dm Dp 

that indicates at which rate two differently filled adjacent boxes equilibrate. The condition r e r M , that reads 
b 3> a/ir 2 exp(— Uo/k B T), is equivalent to the previous one Dm <C D^. The second condition, r e 3> Tb, being 
equivalent to a/b 3> exp(— Uo/k B T)/ir 2 , is automatically satisfied. 

The discrete diffusion equation l|18p is solved by introducing its Green's function: if Wi(t = 0) = Si.o (Kronecker 
symbol), then Wi(t) = Ii(2L>Mt) exp(— 2L>m£); with Dm — Dm /a 2 , can be checked to be the (unique) solution [l9l |. 
In particular, owing to the properties of the modified Bessel functions of the first kind Ii(z), Wi(t) ~ (inDMt)^ 1 ^ 2 
at large t, for any fixed i. This large t behavior will be refined in the next section. 

In the following, we shall also need the box-box, two-time correlation functions of the Wi(t): 

7«(*) = (WitfWjis + t)) - (WitfKWjis + t)). (21) 

Here and in the following section, the averages (•) must be understood as the long-term ones, (-)it- This correlator can 
be computed thanks to the following argument: d/dt Jij(t) = (Wi(s) d/dt Wj(s + t)), and Eq. (TTg|) can be injected 
in this relation, which shows that the jij obey the same diffusion equation as the W, . We recall that (TT*) = 1. Thus 
"fij(t = 0) = Sij, and it ensues that 

jij(t) = Ij-i^Dut) exp(~2D M t). (22) 
This relation will be useful to calculate C(t) below. 
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C. Calculation of the long-term correlation function for a uniform observation profile 

We first consider the correlation function at large times: t > r m . As it will be confirmed in the ensuing subsection 
(devoted to short times, t < r m ), boxes can be considered at equilibrium in the long-term regime because local box 
fluctuations are then decorrelated. The normalized observation profile is I(x) = ^l[-i/2,L/2l- The one-dimensional 
correlation function C(t,L,a) depends obviously on t, but also on the profile width L (it is the core of the present 
work) and on the position a of this profile relatively to the meshgrid (see Figure [5]) . By definition, and after some 
algebraic manipulation using in particular the correlators 7^ , 

C(t,L,a) = [Wi(s)-{Wi)} / 1[-l/2,l/2)(x - a) dx 

i=-oc Jla 

x E [^i( s + *) - ( W i)] / M-W/2](y - a) dy) (23) 

Jja 



] = -oo 

00 00 



= 7 E E I^DMt)e-^l t {a)l 3 { a ) (24) 



L 

l— — OO J — —OC 



with 



r(k+l)a 

l k (a) = 1[-l/2,l/2](x ~ ") dx. (25) 

Jka 

This quantity is in general equal to a or 0, except when k is the index of a box only partially covered by the observation 
profile (see Figure [5]). The prefactor l/L in Eq. (|2"4")l is the product of y/S e s by the square of the normalization of the 
profile, namely (l/L) . 

To go further, we need introducing some notations (Figure From now, all lengths are expressed in units of a 
(thus a = 1 and Dm = Dm)- We set L = i + e with I an integer and e € [0, 1[. The shift a E [0, 1[ is the distance 
of the left-hand-side extremity of the segment to the closest obstacle on its left. We also introduce the integer l m ax, 
indicating the number of boxes fully covered by the observation segment: l max = ^ifa + e< 1 and l max = I — 1 if 
a + e > 1. In addition, two end-boxes are generally only partially covered by this observation segment. Their lengths 
are f3 = 1 — a and 7 = L — l ma x — P- 

I =3 

" max *-* 

— ■— I 1 ! 1 1 1 ► 

a ; ; ; x 



FIG. 5: Definition of the various lengths used in the calculation of Section Till CI The short vertical lines indicate the positions 
of the semi-permeable barriers. All lengths are in units of the mesh parameter a. In this example, the observation segments (in 
red) have length L = 3.6. The number of boxes that it covers entirely is denoted by l ma x- Depending on its position a 6 [0, 1[ 
relatively to the meshgrid, l m ax can be equal to £ = 3 (top) or t — 1 = 2 (bottom). The residues and 7 give the overlap of 
the observation segment with the first and last observed boxes. 

Setting z — 2DMt, and collecting the Bessel functions of same index, the correlation function is re-written as 



(26) 
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Now we can define the average of C(t, L, a) over all possible shifts a: 

C(t,L) = [ C(t,L,a) da 
Jo 

= [ C(t, L, a) da + f C(t, L, a) da. 

JO Jl-s 

In the first (resp. second) integral above, l ma x — £ (resp. lmax = £ — 1). After integration, one gets: 



(27) 



C(t,L) 



(1 - ^)Io(z) + 2^(1 - + (~ + e + e 2 - ^ 



e 3 l e (z) e 3 I t+l (z) 
3 ' L 3 L 



(28) 



Figure [3] shows an example of this long-term approximation (lower curve): it perfectly follows the exact C (upper 
curve) in the region of interest t > r m . Now we use the following large z expansion [l|| where /i = 4j 2 : 



e-'Iiiz) = 



1 



'2-KZ 

1 

/ 2nz 
1 

/ 2nz 



1 - 



/i-1 , (/x- 1)^-9) ( M -1)(m-9)(m-25) 



-fj,/8z 



8z 
1 - 



2!(8^) 2 
2z \4 8z 3 V82 



3!(8z) 3 



1 



/i 



2z V8z/ l\2z 



1 



after re-summing the terms of order 0(1) and 0(1/ z). We defined the polynomial V(X) 
that a; is the adimensional variable of order 1: 



(29) 



X + ±X J . We recall 



/2z 



(30) 



Thus 2z = x 2 L 2 and n/(8z) = j 2 /(L 2 x 2 ). Now at the leading order in 1/z, only the first two terms contribute in 
Eq. (HU) and 



C(x) 



' Ea-^)exp(-fc) 



3=0 

V2 - 1 



2 

2? 



(1 — u) cxp( j)du, 



(31) 



setting u = j/L. Thus C(a;) = r(x) + 0(1/ L) with 



r(x) = erf(l/x) 



- : (l^exp(-l/x 2 )), 



(32) 



as expected for FCS with an uniform observation profile in the case of free diffusion (see Section III Dl) , since at large 
time-scales, in other words at large observation length-scale L, the diffusion recovers a free character. The graph of 
r was displayed in Fig. [TJ The calculation of the next term in the expansion of C in powers of 1/L is derived in 
Appendix [X] It requires in particular to take into account the corrections due to the discrete summation over j. It 
turns out that the terms of order 1/L cancel and that the first non- vanishing term is of order 1/L 2 : 



C(x)=T(x) + ^ + 0(l/L 3 ) 



(33) 



with 



9( x ) = 



-i (1 - exp(-l/x 2 )) + JL exp(-l/.* 2 ) 



(34) 
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Now the decorrelation time till is defined by C(2 V^ lliZi) — l/y/2 and x by r(x ) = 1/V2 for x a ~ 0.520. We recall 

2 

that it means that, at the leading order in L, ti/ 2 = jjy—L 2 (see Section TlIPp . Owing to the work of Wawrezinieck 
et al. [H, we anticipate that 



1/2 



4L> 



M 



-{L 2 -a) + 0(l/L) 



(35) 



where a is the shift in L 2 observed on tj/2 vs L 2 plots in numerical simulations (a is also in units of a 2 ). In other 
words, 



, 1/ ^2V^=,ofl-^) + 0(l/X3 



(36) 



Now 



0(xi/ 2 ) = r(a; 1/2 ) + 



L 2 



0(1/ L 3 ) 



T(x ) + T'(x )(x 1/2 - ar ) + ^ + 0(1/L 3 ). 



(37) 



The second order term vanishes because C(x 1 / 2 ) — r(xo), thus 



2.9 (x ) 
x T'(xq) 



1 



1 2 1 

3^3^ exp(l/^) - 1 
0.49918... 



2xq 



(38) 

(39) 
(40) 



This quantity having been calculated with a high-precision software, it is close to, but definitely different from 1/2. 
However, the value a — 1/2 is an excellent approximation for a square meshgrid and a square observation area. 




2 3 4 

L 2 (units of a 2 ) 



FIG. 6: Decorrelation time, ti/2, in arbitrary unit To, vs the observation domain area, L , for a uniform observation profile 
and Dm — 10 _4 a 2 rJ~ 1 . The continuous line shows the numerical solution of C(ti/2,L) = l/y/2 with C given in Eq. ((28}. The 
dotted line is the result of the asymptotic expansion of ti/2, valid at large L. It crosses the L 2 axis at a ~ 0.5. The numerical 



solution meets the I? axis at V ~ 1.30. 



We have shown in this section that, for a square uniform observation profile, up to corrections of order l/L or 
higher, 



; 1/2 c^IV- 0.499 a 2 ), (41) 
Dm 
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where the length units are again expressed in real units (we recall that a is the mesh parameter). This law is illustrated 
in Figure [5] The abscissa-shift a ~ 0.499 a? is the consequence of the difference between: (i) the correlation function 
C(t), when one takes carefully into account the effects of semi-permeable barriers on diffusion; (ii) and the same 
function T(t) when one only takes into account the long-term effect of these barriers, which results in free diffusion 
with an effective coefficient Dm- This area-shift a amounts to a decrease of the diffusion time t\/2 with respect to 
the zeroth-order approximation t\i<x — 0.0676 L 2 /Dm- The main contribution to a comes from the term 1/3 in the 
brackets of Eq. ([39|) . the second term being approximately equal to 0.06. By inspection of Appendix [A| this term 1/3 
has 3 different contributions: 1/3 = 4(— 1/12 — 1/6 + 1/3). The first two negative contributions come respectively from 
the short-time corrections to the approximation of the Bessel functions (materialized by the polynomial V, Eq. (|A5|) ) 
and the approximation of the sum on j by an integral (Eq. (|A6[) ) . The only positive contribution to a, that gives it 
its final positive sign of interest here, come from the two extremities of the observation segment of length L (segments 
f length (3 and 7 in Fig. [5J) : there, the local in-box correlations fall off much more rapidly than the global correlations 
between boxes, which reduces the overall decorrelation time ii/2- Beyond er, the next-order corrections to t-1/2 can 
in principle be calculated by expanding the above calculations at order 1/L 3 . It would allow to estimate for which 
values of L 2 the approximation (|4ip is valid. However, a rapid inspection of Fig. [5J shows that the approximation is 
excellent whenever L 2 > 6 a 2 or L > 2.5 a. At shorter L 2 , the next-order corrections cannot be neglected and must 
be taken into account when exploiting experimental data. This issue will be addressed in the Discussion (Section HV|) . 
together with the additional effect of the short-term decorrelation neglected so far. 



D. Calculation of the short-term correlation function 



Our purpose here is to understand the short-time (t < r m ) behavior of the correlation functions, in particular the 
first fall-off of C(t) observed in Figure [31 In contrast to the previous calculations, the short-term averages (-) s t must 
be taken into account here. The interest of working with a uniform observation profile is that only end-boxes, of 
lengths (3 and 7, will contribute at short times, which renders easier the understanding of underlying physics. 

First of all, we consider the box of index i, the boundaries of which are approximated by perfectly reflecting ones. 
This approximation is correct at times t < r m -C r e because exchange of matter between adjacent boxes is negligible. 
This box is observed with a uniform observation profile of intensity 1/L and length 77, one extremity of which coincides 
with a box boundary, situated without loss of generality at x = 0. This case includes both the two end-boxes (with 
77 = [3 or r) — 7) and internal boxes that are entirely covered by the observation profile (with rj — a). This observation 
profile is written l(x) = l[o, v ](x) / L. The total box filling, Wj, is constant, but can be different from 1 here since it 
fluctuates slowly at large times. The particle density at the position x in the box, p(x,t), is written ^ L ip i (x,t), this 
last function ipi taking into account short-time local fluctuations in the box. It is a Poisson random variable of mean 
(^i)st = 1 and mean square (ipi(x,t)tlji(y,t)) st = 1 + ^S(x - y). 

Like free diffusion (Section III D|) . this case is tackled by Fourier analysis, which requires the calculation of the 
Fourier modes of the diffusion equation in this special geometry. We shall not detail these calculations here, which 
can be found in text-books (20l |. but simply give the final result. The short-term contribution of this box i to C(t, L, a) 
is 



LW 2 



dx 



dy i{x)i{y) [(ipi(x,s)ipi(y,s + t)) at - (ipi) a t(ipi) at] 



(42) 



where 



G v (t) = - 



1 1 00 1 

£+4 y - 



sm 



V 

nir— 

a 



exp 



D 2 2, 
— £-71 n • t 
a 2 - 



(43) 



The prefactor L = %/Srff m Eq. (j42|) comes again from the overall normalization of C(t, L, a). First of all, as anticipated 
in Eq. (fT^jl . these correlations fall off exponentially for times larger than 77 = a 2 /(7r 2 I? At ); G v (t) decreases from 1 at 
t = to rj/a when f > 77. If 77 = a, then G v (t) = 1 is constant. Internal boxes do not participate to the short-term 
decorrelation. In addition, consider two different boxes, with indices i and j, i ^ j. Since inter-box currents vanish 
at the time-scale under interest, the short-term correlations, (ipi(x, s)ipj{y, s + t)) st — 1, are also negligible. Therefore 
the only contribution to short-term decorrelation comes from the auto-correlation of both end-boxes. 

Now we must take the long-term average (-}it, which amounts to replacing Wi by (Wi)u — 1. We set again a = 1. 
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The short-term contribution of both end-boxes to C(t) is then 



L 



The exact shape of this contribution appears in Figure [3l at times t < r m (after integration over a). One can remark, 
from comparison of the algebraic expressions pi?)) - where e~ z Io(z) ~ 1 since z — 2E>Mt — 2t/r e « 1 - and (|44[) . 
or from direct inspection of the figure, that our short- and long-term correlations functions coincide at intermediate 
times t rs -' Tin . This is £lh additional validation of our approach. 



p + i 

L 

/3 2 + 7 2 



if t = 



j if t>r M 



(44) 



E. Calculation of the long-term correlation function for a Gaussian observation profile 

Here we focus again on the long-term correlation function, but in the Gaussian case. The short-term one can be 
treated in the same way as in the square case, with the additional complication that all boxes will contribute to it. 
The final result will be that the short-term contribution will also vanish exponentially for times larger than r M , since 
this is a property associated with the boxes of size a and not with the precise shape of the observation profile. By 
contrast, the numerical pre- factors of the decaying exponentials will be different from those of Eq. ([43]) (the squared 
sines). However our main goal in the present paper is to tackle the calculation of the shift a, which only is of practical 
interest @. This requires the sole long-term correlation function C. Its calculation follows the same route as in the 
previous square case, in Section UlI CI We shall see that the dominant contribution comes now from the approximation 
of the discrete sum over boxes by an integral. The long-term diffusion is enhanced by the rapid diffusion in boxes 
that rapidly equilibrates the local particle densities. 

The normalized one-dimensional Gaussian illumination intensity profile is now (see Section III Bp : 



I(x) 



V2 



■ exp 



2x^_ 

Hi 2 



(45) 



where the waist, w, plays the same role as the length L in the previous uniform case. Again, this profile can be 
centered anywhere relatively to the meshgrid. The function C(t,w,a) is defined as in the square case (Eq. (f24l) ): 



oo oo 



C(t, w, a) 



J2 E I^D M t)e- 2D "%(a)l 3 (a) 



(46) 



l= — OO J = — oo 



with 



1k{a) 



exp I g i <> ; • 



(47) 



The prefactor is the product of y/S c s by the square of the normalization of the Gaussian, ^ ^F^ . We first 
integrate on the position of the center a of the intensity profile: 



C(t,w)= / C(t,w,a)da=^- V V 7 i _ i (2 J D M t) e - 2Z3M */C ii , 
Jo ■ ■ 



%=■ — OO J = — OO 



(48) 



where 



Xi(a)Xj(a) da 



dx dy / da exp 
Jj Jo 



:[{x - af + {y - af 



(49) 



Now we develop the exponential in powers of 1/w and a: (x — a) 2 + (y — a) 2 = h(x — y) 2 + 2[a — h(x — y)] 2 , thus 



exp ( ~[(x - a) 2 + (y- a) 2 }) = exp ( -- L [(x - y) 2 + (x + y) 2 } 



+ 0(l/w 3 ). 



(50) 
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Here and in the sequel of this section, one must keep in mind that the calculation will end with a summation over 
i and j (Eq. (|46p ). Therefore, the term in the square bracket proportional to (x + y) vanishes when summing over 
i and j, by parity. In addition, the exponentials in Eq. f50|) are Gaussian weights of the order of exp(— (i + j) 2 /w 2 ) 
or exp(— (i — j) 2 /w 2 ). Thus, when summing over i and j, any factor of order (i + j) n or [i — j) n will have a final 
contribution of order w n that must be taken into account when developing any expression in powers of 1/w. This 
is precisely the case for the factor [x + y) 2 in the previous expression, which will eventually give a contribution of 
order w 2 because x « i and y w j in Eq. (|49p . For this reason, the third term in the bracket, proportional to a 2 , 
also vanishes at order 1/w 2 when summing over i and j. Consequently, at order 1/w 2 of interest here, there is no 
dependency of C(t,w) on the observation profile position, a. The integration over a to compute ICp is trivial. This 
remark, together with the forthcoming calculations, will confirm the numerical observation in Ref. p (Fig. 10B): at 
large w, decorrelation times t 1 / 2 become independent of a. 

The Bessel functions are treated again using Eq. (|29p and at order 1/w 2 , still using z — 2D^t, 



C(t,w) = ^z~7t,^1 ex P 



TT W\ Z 



u - s 

2z 



1 ij- i) 
1 -| -p y A L 

2z \ 2z 



h3 

I dx j dy exp^- — [(x-y) 2 + (x + y) 



(51) 



At the leading order in powers of 1/w, defining again x = y2z/w, we get C(x) = T(x) + 0(l/w) with 
V2 1 



T(x) = 



7T Wy/ Z 
1 



dx / dy exp 



TTXW^ 



— oo •/ — oo 
oo 



l + x z 



du I dv exp ( — j ( y2 + y2 ) 



1/2 



2z 

u 2 1 

2z w 



[u = y — x; v = y + x] 



(52) 



by Gaussian integration, as expected for FCS with a Gaussian observation profile (see Section l"IID|) . In this case, 
xo = 1. 

The corrections of order 1/w 2 to T(x) come from: (i) the term that is explicitly in l/(2z) = l/{x 2 w 2 ); (ii) the 
discrete summation over i and j, as in Eq. (|A6[) . The case (i) is easily handled by Gaussian integrals. It gives a 



correction 



V2 1 
32 up 



at order 1/w 2 . The case (ii) requires a finer analytical investigation to be sure to track all the terms 



of the relevant orders. It is treated in Appendix [51 and gives another correction — -ori- One finally gets 



C(x) = T{x) + 4, + 0{l/w z ), 



with g now independent of x: 



Finally, as in the previous uniform case, 



V2 
32 

2.9 



ur 

V2 
24 



V2 
96 ' 



— ~ 0.0833. 



(53) 



(54) 



(55) 



x T'(x ) 

This quantity can now be compared to the results obtained by numerical simulations by Wawrezinieck et al. 0], by 
relating it to the quantity X 2 computed there: X 2 = na ~ 0.262 (in their case, the meshgrid parameter is r = a/2). 
Owing to the statistical errors inherent to any numerical method, this quantity had been over-estimated in previous 
publications and had been found to be larger than 1. However, the present analytical calculations confirm that 
the phenomenon pointed out by these numerical investigations enables the detection by FCS of underlying barrier 
meshgrids impeding long-term diffusion. We have also computed numerically the exact tin vs L 2 plots, as shown in 
Fig. El by solving the equation C r (t 1 / 2 , L) — l/y/2, where C{t, L) is given in Eq. (|4"5|). They resemble those of Fig. [6] 
One notices that in this case, the order-two approximation 

1 



&1/2 



AD 



-(w 2 -a 2 /12) 



(56) 



is excellent as soon as w 2 > 2 a 2 . Below this limit, t\i% is sensibly smaller than its affine approximation and vanishes 
at L 2 ~ 0.362. 
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FIG. 7: Decorrelation time, ti/2, in arbitrary unit to, vs the observation domain area, w 2 , for a Gaussian profile of waist w 
and Dm = 10 _4 a 2 r ( j~ 1 . The continuous line shows the numerical solution of <5(i]/ 2 ,w) = l/v2. The dotted line is the result 
of the asymptotic expansion of ti/2, valid at large w. It crosses the w 2 axis at a = 1/12 ~ 0.0833 (see text). The numerical 
solution meets the w 2 axis at w 2 ~ 0.362. 

IV. DISCUSSION 

In this work, we have confirmed, by exact analytical arguments, previous numerical findings of P.-F. Lenne and his 
collaborators [5, 6]: FCS at variable observation area gives a typical signature of two-dimensional diffusion hindered 
by a meshgrid of semi-permeable barriers. The decorrelation time t\/2 is close to an affine function of the observation 
area w 2 and its experimental determination allows that of the sub-wavelength mesh parameter. This technique has 
indeed been shown in the same references to be applicable to the characterization of the diffusion properties of proteins 
embedded in the plasma membrane of live cells. In this case, the meshgrid can be constituted of cortical filaments of 
the cytoskeleton interacting with the cytoplasmic protuberance of proteins [H, [l3| . 

For the sake of pedagogy, we have first presented our calculations in the case of square, uniform, observation profiles 
before focussing on Gaussian ones. Conclusions are qualitatively similar in both cases. We recall that w and L/2 
(or w 2 and L 2 /4) played essentially the same numerical role in the case of free diffusion. To what respect does this 
similarity survives on a meshgrid? First of all, by visual inspection, one can see that Figs. [6] and [7] are quite similar. 
Being equivalent to free diffusion, their long-time behavior is indeed identical, providing that one identifies w 2 and 
L 2 /4. Besides, at shorter times, the ratio of cr's is 0.499/0.083 = 6 «4, and the ratio of points where the numerically 
calculated ti/2 vanish is 1.30/0.362 = 3.59 w 4. Therefore the correspondence L/2 <-> w is essentially preserved when 
going from the uniform profile to the Gaussian one. To this respect, there is another interest in having studied the 
square uniform case: circular, uniform profiles cannot been tackled within our formalism. But we expect them to be 
in-between the Gaussian ones (which have the circular symmetry), and the square ones (which possess the uniform 
character), and thus to resemble the two previous ones. 

Now, we restrict the discussion to Gaussian profiles of experimental interest. We have proven that the decorrelation 
time £1/2 (at which the two-dimensional correlation function gets reduced by a factor 2), displays an affine behavior, 

h /2 = -±-(w 2 -a 2 /12), (57) 

provided that the waist, w, and the mesh parameter, a, satisfy w 2 > 2 a 2 . As compared to the numerical simulations 
previously mentioned, this work provides an accurate value of the parameter a equal to 1/12 in this equation. The 
numerically calculated value was over-estimated by a factor larger than 4, thus a priori leading to an under-estimate 
of the mesh parameter a by a factor 2. 

The former over-estimate of this parameter 1/12 does not go into the right direction for a good experimental 
extraction of the mesh size a. Indeed @, @|, to have a good estimate of the shift a, one has to go down to observation 
length-scales only slightly above y/a ~ a/3. When working at laser wavelength A, this imp lies a/3 > A/2, i.e. 
a > 600 nm, while in general mesh sizes are supposedly not larger than 100 nm 0, HH, l2l| . This problem can 
be bypassed by employing optical techniques overcoming the Rayleigh limit, such as nanometric apertures [f| or 
Stimulated Emission Depletion (STED) [22|, both of which a priori enable one to go down to deci- nanometric length- 
scales. In the latter case, the profile is also close to Gaussian .22|. But another complication arises when working 
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at this scale: near a, t\/2 is not an affine function of w 2 , as shown in Figure. [7j Nevertheless, it can reasonably be 
approximated by an affine function that now crosses the w 2 axis at w 2 ~ 0.362 (corresponding to X 2 = it 0.362 ~ 1.14, 
with the notations of Refs. 0, It is this value of X 2 that must be used for these small values of L 2 . This last 
argument explains the value X c sa 1 found numerically in these references, where £1/2 is fitted linearly without taking 
into account the curvature of the numerical curve observed in Figure. [71 finally re-conciliating numerical and analytical 
approaches. 

We should also emphasize the following important points. Our calculations rely in great part on the effect of the 
first fall-off of C(t) due to microscopic decorrelation, around r M . The calculation of £i/ 2 takes into account this fall-off 
because it reduces the value of £1/2 with respect to free diffusion. But to see this effect on experimental data, one 
must be sure that the FCS curve C(£) is acquired on a time range, [t m i n ,t max ], that does include this fall-off. In 
other words, to calculate when C(£) is equal to 1/2, one must know the reference values C(£) = at large times and 
C(t) — 1 at short times. Thus one must choose t m in strictly lower than £ M . Diffusion coefficients are smaller than 
0.1 /im 2 s -1 for membrane proteins in live cells, and a is larger than 30 nm, as discussed just above. Thus is much 
larger than 1 ms. A value of t m in of 0.1 ms should be sufficient to ensure that short-term decorrelations are correctly 
taken into account in all live-cell plasma membranes. 

Furthermore, fluorophorc photo-physics are responsible for additional noise in FCS profiles at very short times. The 
most common cause for flickering in the fluorescence intensity is the transition of dyes to their first excited triplet 
state, where they stay a relatively long time before relaxing to their ground state [23| . If the contribution of this effect 
is superimposed to the short-time decorrelation addressed in the present study, the first fall-off will be over-estimated 
and the corrections to £1/2 will contain a systematic error. However, this noise decorrelates after a few fis [HI, and it 
will not interfere with the mechanisms of interest here if one ms. But this remark shows that one 

must be sure that all sources of noise specific to fluorophores are completely decorellated at £ m m in order to analyze 
correctly the diffusion on a meshgrid. This remark should help experimentalists to orient their fluorophore choice. 

In the previous section, we have also noticed that as far as the Gaussian observation profile is concerned, the 
dependency of £i/ 2 on the observation profile position a vanishes. At the order of 1/w considered here, which is 
adapted to study the large w behavior of t 1 / 2 , <J does not depend on a. This point had already been observed in 
numerical simulations [5, 6]. However, if we had worked at a higher order in 1/w, we would have found that a has 
corrections of order 1/w or higher that do depend on a. Indeed, at small w, it is clear on the simulation results that 
the plots of t 1 / 2 vs w 2 depend on a. By contrast, in the case of a uniform profile, if a is fixed, then £i/ 2 depends on 
a (calculations not shown) because g a {x) depends explicitly on a. We have chosen in our calculation to average over 
a before calculating a: g(x) = L g a (x) da. At the experimental level, this corresponds to averaging FCS correlation 
profiles over realizations (at fixed w) before extracting £i/ 2 . However, we could have chosen instead to calculate 
£1/2(0;) (or equivalently o~ a ) for each value of a before averaging £1/2(0) over the different values of a. This would 
mean experimentally that one calculates £1/2 for each realization before averaging over realizations. But a a is linear 

in g a (xo) (Eq. 1|38[)) and L o~ a da — 2 f g a (x) da / (xqT' {xq)) = o. Thus the overall result does not depend on the 
way of averaging over realizations. 

Another issue that we have not debated so far is the contribution of the confined microscopic diffusion in boxes. 
If Dm/D^ is very small, then r M ~ a 2 /D^ is also very small and the first fall-off in Figure [3] is pushed away to the 
left. It will have a significant influence on £i/ 2 only when L 2 is very close to a. One then expects its contribution to 
be of order L 2 /D^. Since Dm/D^ <C 1, this contribution is completely negligible as compared to the one previously 
discussed. In contrast, if Dm < D^, but their ratio is not vanishingly small (typically Dm/D^ ~ 1/10), one can still 
speak of diffusion on a meshgrid. But the short-term contribution cannot be neglected anymore in the calculation 
of £1/2, as it is visible in numerical simulations 0, @. This makes the estimation of a 2 more difficult, because the 
point a where the long-term contribution of £i/ 2 crosses the L 2 axis is now hidden in the short-term contribution 
(see Refs. [1, Q). Even though it is outside the scope of the present paper to tackle this issue, the calculation of 
this effect is feasible in principle, starting from the full correlation function. In addition, the microscopic diffusion 
constant, D M , is also theoretically accessible from experiments, because is a priori measurable at the first fall-off of 
the correlation function when L is close to a (even if the FCS signal appears to be quite noisy for these time values). 
Once is known, if it appears to be large as compared to the measured value of Dm, then one can forget this issue. 
In the converse case, the knowledge of both D^ and its theoretical contribution to £i/ 2 enables one to subtract this 
contribution from £]/ 2 and thus to extract a good estimate of a. This work remains to be done. 

We must also discuss the fact that at the experimental level, the underlying meshgrid has no reason to be a periodic, 
square lattice. It is more likely to look like a random meshgrid with cells of variable sizes and shapes. One cannot 
expect from an analytical approach to solve such a complex issue. Only numerical calculations can provide reliable 
results. However, it is reasonable to expect that the mean mesh diameter will be very close to the size a measured via 
the procedure discussed here. Furthermore, the present work can be generalized to three-dimensional diffusion in a 
similar meshgrid of obstacles defining three-dimensional, semi-permeable cages (e.g. diffusion in a gel-like cytoskeletal 
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mesh). One has simply to replace Cjt T ){t) = C 2 (t) by C/( r )(t) = C 3 (t) in our analysis. The diffusion time, ti/ 2 , is 
then defined by Cji T ){ti/2) = 1/2 v2- All the argumentation remains valid. 

Finally, the present technique is not able to tackle isolated domains where membrane proteins are temporarily 
trapped and diffuse more slowly, as described in Refs. @, @]. It has been shown numerically by P.F. Lenne and his 
collaborators that in this case, t\/2 still displays an affine behavior, but that a becomes negative. The reason for the 
failure of our approach is that the potential U(x,y) cannot be written as a sum of one-dimensional potentials in this 
case. What happens if we write U(x,y) — U\(x) + U±(y), where the U\{x) are one-dimensional, periodic potentials, 
alternating zones of low and high potentials and of slow and rapid diffusion coefficients, D^l Then diffusers encounter 
bands, both in the x and y directions, where the potential is higher (or lower) and the diffusion slower. At the 
intersections of these bands, there are square regions where U is even higher (or lower) and where diffusion is even 
slower. These isolated, square regions could mimic those perviously mentioned and this problem could be tackled 
from our point of view. However, the quasi-static box approximation is not valid anymore in this case because we 
want b to be of the same order of magnitude as a, thus making unrealistic the condition b -C a. Here also, work 
remains to be done in the future. 
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APPENDIX A: EXPANSION OF C IN POWERS OF 1/L (UNIFORM CASE) 



In this appendix, we identify and calculate the different corrections of order 1/L or 1/L 2 to Eq. (|3"2")) . First of all, 
this expression (f3"2")l comes from the approximation of a discrete sum, the first two terms in Eq. (|28p. by an integral: 
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where f(v) = (1 — v) exp(— v 2 /x 2 ) ^1 + L l xi V ^^H^. The first term B in Eq. (|A2|) is the approximation by the 
trapeze integration method of the integral: 
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and the residue, R, is given by the standard error of the trapeze method on each interval \j/L, (j + 1)/L]: 
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Finally, the three remaining terms in the brackets of Eq. (j28|) of order 1/L contribute to C(x) at the order l/L 2 



IHxL 2 



--(1 - exp(-l/x 2 )) + (e + e 2 ) exp(~l/x 2 ) 



+ <%)■ 



(A7) 



Thence 



C(x) = Y(x) + 



^/ttxL 2 



^(l-exp(-l/x 2 )) + ^exp(-l/x 2 ) 
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APPENDIX B: CORRECTIONS OF ORDER 1/w' 2 DUE TO THE DISCRETE SUMMATION IN EQ. (^Tj) 

(GAUSSIAN CASE) 

In this appendix, we calculate the corrections of order 1/w and 1/w 2 to C(x) = in the Gaussian case, due to 

the discrete sum in 



D(t,w) = Z — T E ex P 
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dx / dy exp ( -[{x — y) 2 + (x + y) 2 ] 
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This expression is that of C(t, w) (Eq. (|5ip ) where all terms being explicitly of order 1/w 2 or 1/z have been omitted. 
However, the summation cannot be replaced by an integral without taking care of the ensuing errors because they 
will also appear to be of order 1/w 2 . We set again the new variables u = y — x and v = y + x, and 
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The 1/2 prefactor comes from the Jacobian of the substitution of variables. First we expand the integral over u in 
powers of 1/w: 
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The next step consists now of integrating over v while keeping track of the relevant terms in the expansions in powers 
of 1/w. Tedious but straightforward calculations lead to 



D(t, W ) = — 7=yz ex p - » ) ex p 



2z 



1 



[(»"j) 2 + (^ + .?) 2 



(i+j + l) s 
3w 2 



6txi 2 6w 2 ^ 



(3 -if 



0(l/w 3 



(B4) 



Following the same route, we expand 



E(t,w) 



V2 1 



dx / dy exp 



(x- yf 
2z 



exp o[[x-y) + (x + y) ] 



(B5) 



We get an expression comparable to D(t, w), except that in the last term of the square bracket of Eq. (|B4I) . both w 2 
must be replaced by a 2 , where a 2 = + is of the order of w 2 . Now E(t, w) is a mere Gaussian integral that 

is equal to T(x). By subtracting from Eq. (|B4[) its equivalent for E(t, w), we finally get: 



D(x,t)=T(x)-^^ + 0(l/w 3 ). 



(B6) 
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These tricky calculations, as well as those of Appendix [5] have been checked using a numerical integration software. 
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